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PARAMETRIC STUDY OF RELAMTNARIZATION OF 
TURBULENT BOUNDARY LAYERS ON NOZZIE WALLS 

J. P. Kreskovsky, S. J. Shamroth, and H. McDonald 
United Aircraft Research Laboratories 


SUMMARY 


By means of an extensive con5)arison between theoretical predictions and 
experimental data this investigation asRPggeg the accuracy of a boundary layer 
procedure to predict the effect of large streamwise accelerations upon 
initially turbulent boundary layers. The boundary layer procedure under con- 
sideration is based upon simultaneous solution of the boundary layer partial 
differential equations and the integral turbulence kinetic energy equation and 
has previously been shown to yield satisfactory predictions for a leirge number 
of boundary layers in the subsonic to low hypersonic Mach number range. The 
results of the present investigation show the ability of the procedure to 
accurately predict properties of boundary layers subjected to large streamwise 
accelerations. Hie procedure is used to conduct a parametric study of the 
effect of free -stream turbulence, heat transfer, Reynolds number, acceleration, 
and Mach number on boundary layers in supersonic nozzles to assist in the 
design of a quiet tunnel. Results of the investigation show that, even in the 
presence of moderate free-streeun turbulence levels, the boundary layer in the 
approach section of the quiet tunnel nozzle re laminar izes and becomes thin 
enough to be removed by a small slot in the nozzle wall. Furthermore, the 
calculations indicate that it should be possible to maintain a laminar boundary 
layer for the entire length of the supersonic portion of the quiet tunnel 
nozzle. 


INTRODUCTION 


The location of boundary layer transition is an important consideration 
in the successful design and operation of supersonic and hypersonic vehicles. 

In general, it is expected that the maximum wall heating at any point in the 
flight path will occur in the region where the boundary layer undergoes 
transition from a laminar to a turbulent state and the insurance of vehicle 
structural integrity requires an accurate prediction of both the transition 
location and the peak wall heating rate either by theoretical or experimental 
means. In addition to heat transfer considerations, boundary layer transition 
may play a significant role in determining the skin friction drag and vehicle 
wake flow field. At present the only method of determining the transition 
region for the complex three-dimensional flow fields developing on hypersonic 
vehicles other than flight testing is through accurate wind tunnel testing. If 
a model is tested in the proper environment, it should be possible to obtain a 
flow pattern in the wind tunnel corresponding to the flow pattern which would 
occur under free flight conditions. Two obvious free flight quantities which 
must be matched in the wind tunnel test are Mach number and Reynolds number . 

A third quantity which must be matched particularly in regard to determining 
the transition region is the free-stream disturbance level. It has long been 
recognized that the free-stream turbulence level has a dominant role in 
determining transition (e.g., refs. 1 and 2) and, therefore, it is important 
that the wind tunnel tests be characterized by the free-stream disturbance 
level expected to be present under flight conditions. 

The problem of obtaining a specified free-stream disturbance level in a 
high Mach number wind tunnel is one which requires careful analysis during the 
tunnel design period as indicated by the fact that transition measurements 
carried out in different wind tunnel facilities have been inconclusive and even 
contradictory. For example, Pate and Schueler (ref. 3 ) Pate (ref. h) have 
shown that correlations based upon transition data from one supersonic wind 
tunnel do not necessarily correspond to correlations based upon data from a 
different tunnel. Upon examination of the data, Pate and Schueler (ref. 3 ) 
showed that for a large amount of data a direct relationship exists between the 
location of transition and the free-stream aerodynamic noise emanating from 
turbulent tunnel wall boundary layers, thus indicating that at high I'fech 
numbers the acoustic mode is the dominant free-stream disturbance mode and that 
the source of the acoustic disturbance is the transitional and turbulent 
portions of the tunnel wall boundary layer. Bertram and Beckwith (ref. 5 ) have 
also demonstrated that a considerable portion of the noise found in a supersonic 
tunnel emanates from the turbulent tunnel wall boundary layer. Thus, the 
development of a quiet tunnel more representative of free flight conditions 
would require a laminar tunnel wall boundary layer in all regions which could 
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effect the test section of the tunnel. One method of suppressing the turbulent 
boundary layer is to cause the normally turbulent boundary layer which develops 
in the nozzle inlet section to revert to a laminar-like state by means of high 
strearawise acceleration (ref. 6). Furthermore, wall cooling may delay 
transition in the nozzle approach section and may suppress transition of the 
laminarized boundary layer in the supersonic section of the nozzle once the 
strearawise acceleration is removed. Reynolds number and Mach number affects 
are also expected to influence the re laminar izing boundary layer. Thus, an 
analytical procedure which could predict the effects of high strearawise 
acceleration, wall cooling, Reynolds number, and Mach number, as well as the 
effect of free -stream turbulence on transition and relaminarization, would be 
a tool of considerable assistance in the design of the quiet wind tunnel. 

Presently, there exist three possible approaches for developing a 
boundary layer theory to predict traiisition and relauiinarlzation : the semi- 

empirical approach, the stability theory approach, and the turbulence kinetic 
energy approach. Typical examples of the semiempirical approach for predicting 
transition location are discussed by Hairston (ref. 7) and such approaches 
could also be used for predictions of relaminarization; however, semiempirical 
procedures cannot be used with confidence to predict transition for flow 
conditions significantly different from the flow conditions of the correlating 
data. In addition to the limitations imposed by data correlations, most 
semien 5 >irical theories are limited further by the crude assumption of 
instantaneous transition. Although the instantaneous transition assumption may 
be acceptable if flow conditions in the vicinity of transition are not of 
interest, the assumption clearly is unacceptable if predictions in the vicinity 
of transition are required. The instantaneous transition assumption has been 
relieved in the work of Harris (ref. 8) who assumes a finite length transition 
region in which the total shear stress is a linear combination of the laminar 
shear stress and the fully-turbulent sheea* stress. This model has successfully 
predicted the development of transiticaial boundary layers in the low hypersonic 
Mach number regime and represents a significant improvement over instantaneous 
transition models. This model could be extended to relaminarization, however, 
since it heavily depends upon empiricism and was not intended to predict the 
effect of phenomena such as free -stream turbulence, Mach number, wall tempera- 
ture, etc. upon transitional behavior, it does not serve as a general boundary 
layer transition or relaminarization prediction procedure. The second possible 
approach which might lead to an analytical prediction of transitional boundary 
layers is based upon stability theory in which the flow is divided into a mean 
flow whose stability is the subject of the investigation and a superimposed 
disturbance of specified form. Attempts to predict the initiation of 
transition from classical stability theory (for example, ref. 9) have had only 
qualitative success and this approach has yet to be applied successfully to 
the problem of transitional boundary layer development. Furthermore, the 
application of stability theory to relaminarizing boundary layers is not obvious. 
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The third approach for predicting the behavior of transitional boundary- 
layers is based upon the solution of the turbulence kinetic energy equation. 

The turbulence kinetic energy equation has proven to be a useful tool in 
predicting the behavior of a wide variety of turbulent boundary layers 
(e.g., refs. 10 and ll) and has been extended into the transitional regime by 
Glushko (ref. 12), Donaldson, Sullivan, and Yates (ref. 13), and McDonald and 
Fish (ref. lU). At the present time, neither the work of Glushko (ref. 12) or 
Donaldson, et al. (ref. 13) has been developed into a practical prediction 
procedure. However, the approach used by McDonald and Fish (ref. l4) has 
proven capable of accurately predicting the beha-vior of a wide variety of 
transitional boundary layers in the subsonic and low supersonic Mach number 
regimes and this approach has been extended by Shamroth and McDonald (ref. 15) 
to the low hypersonic Mach number regime. 

In the present report the ability of the McDonald -Fish procedure (ref. lU) 
to predict the effect of strong favorable pressure gradients upon an initially 
turbulent boundary layer is assessed and the procedure is then used in 
conjunction with various inviscid flow calculation procedures to assist in 
developing a design for the NASA pilot model quiet tunnel (ref. l6). In the 
quiet tunnel the nozzle is preceded by a constant area entrance duct. The 
entrance duct is long enough so that a turbulent wall boundary layer is present 
at the junction between the entrance duct and the nozzle itself. In the sub- 
sonic portion of the nozzle the flow is subjected to a strong favorable pressure 
gradient which is expected to cause the initially turbulent boundary layer to 
relaminsurize and become very thin. Slightly upstream of the throat the thin 
boundary layer is removed by a suction slot. A new boundary layer is formed 
at the suction slot which grows along the tunnel wall downstream of the throat 
region. This new boundary layer is initially laminar and the favorable 
pressure gradient inhibits transition. Eventually at some station downstream 
of the favorable pressure gradient the boundary layer may undergo transition 
and become turbulent. 

In the present report a series of comparisons between theoretical 
predictions and experimental data for initially turbulent boundary layers under 
the influence of strong favorable pressure gradients are presented, "niese 
comparisons demonstrate the ability of the McDonald -Fish procedure to predict 
the effect of strong favorable pressure gradients on an initially turbulent 
boundary layer. The theoretical procedure is then used to calculate wall 
boundary layers for the pilot model quiet tunnel. In addition to these boundary 
layer calculations, results from inviscid flow predictions for the pilot tunnel 
geometry are presented and compeared to predictions made independently by NASA. 
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LIST OF SYMBOLS 


Cf 

F 

G’ 

H 

H* 

K 

k 

L 


Lref 

I 

■^00 

M 

Pr 

PTiji 

P 

Q 


structural coefficients of turbulence 
specific heat 
skin friction coefficient 
sublayer damping factor 
dimensionless stream function 
dimensionless teni5)erature ratio 
shape factor 

incompressible shape factor 
acceleration parameter 
thermal conductivity 
dissipation length 
reference length 
mixing length 

wake value of mixing length 

Mach number 

Prandtl number 

turbulent Prandtl number 

pressure 

heat flux 

turbulence kinetic energy 



radius 


turbulence Reynolds number 

layer averaged turbulence Reynolds number 

Reynolds number based upon momentum thickness 

Stanton number 

static temperature 

total temperature 

streamwise velocity 

friction velocity 

transverse velocity 

cross flow velocity 

streamwise coordinate 

transverse coordinate 

dimensionless transverse coordinate 

indicator equal to one for axisymmetric flow, 
zero for two-dimensional flow 

density ratio 

intermittency factor 

ratio of specific heats 

boundary layer thickness 

sublayer thickness 

displacement thickness 


reference length 


€ 


turbulence dissipation 


T] dimensionless transverse coordinate 

9 momentum thickness 

U viscosity 

V kinematic viscosity 

v,p kinematic eddy viscosity 

p density 

T shear stress 

X dimensionless longitudinal distance 

01 >02 >03 integral functions (see Eqs. (A-I9) through (A -21)) 

Subscripts 

e boundary layer edge condition 

00 free-stream condition 

w wall condition 

t throat 

Superscripts 

— average quantity 

' fluctuating quantity 



THEORY 


The UARL prediction procedure is based upon a simultaneous solution of the 
boundary layer partial differential equations of conservation of mass, momentum, 
and energy. The partial differential equations are used in conjunction with 
an integral turbulence kinetic energy equation and a turbulence structure model 
to pr’edict boundary layer development. The procedure has been used to predict 
a wide variety of boundary layers with Mach numbers varying from the incompress- 
ible to the low hypersonic regime. The procedure can include the effects of 
heat transfer and wall transpiration and, in addition, has a multiple species 
and chemical reaction capability. In brief, for a single component flow the 
procedure solves the governing momentum and energy equations by a Hartree- 
Womersley approach in which streamwise derivatives are replaced by finite - 
differences, the coordinate normal to the wall is nondimensionalized and a 
stream function is introduced. The resulting momentum and energy equations 
are third order and second order nonlinear ordinary differential equations, 
respectively, in the transverse coordinate. At each streamwise station the 
nonlinear coefficients are estimated from solutions at the previous stations 
and the resulting linearized equations are solved as two point boundary value 
problems. The resulting solutions are used to obtain better estt mates of the 
nonlinear coefficients and the entire process is repeated until two solutions 
agree within a specified tolerance. When two successive solutions agree, the 
equations are considered solved at that streamwise station and the solution 
proceeds to the next streamwise station. 

The integral turbulence kinetic energy equation is coupled iteratively 
to the momentum and energy conservation equations. At specified points in the 
calculation which iterate upon the nonlinear coefficients of the momentum and 
energy equations, the integral turbulence kinetic energy equation is solved to 
determine the turbulent contribution to the transport coefficients. The 
integral turbulence kinetic energy equation is solved even during laminar flow 
and in these cases the predicted turbulent transport properties are negligible 
compared to the laminar transport properties. Gradually as the boundary layer 
thickness increases the predicted turbulent transport becomes comparable and 
then much larger than the laminar transport (except, of course, in the sublayer 
region) indicating transitional and fully-turbulent flow. Similarly in the 
case of a fully-turbulent flow subjected to a strong favorable pressure 
gradient, the predicted turbulent contribution to the total transport may 
decrease leading to relaminarization and development of an effectively laminar 
boundary layer. It Is Important to note that in the UARL procedure transition 
and relaminarization are natural occurrences of the turbulence model rather 
than being triggered by some semiempirical criteria such as momentum thickness 
Reynolds number. In its present form the procedure does not include any effect 
of Taylor -Goertler vortices which if present may induce transition. However, 


8 



if information were available which could relate two-dimensional mean flow 
conditions to the formation of the vortices and to the subsequent development 
of the Reynolds stress tensor, it seems reasonable to assume that this infor- 
mation could be incorporated within the framework of the turbulence kinetic 
energy equation and models such as the present, which are based on energy 
considerations, could be extended to include the effect of the vortices. 

A detailed description of the governing equations and the turbulence model 
has been presented in refs. l4 and 15- The description of the method presented 
in ref. 15 is repeated in Appendix A of the present report. 
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COMPARISON WITH EXPERIMENTAL DATA 


The primary purpose of the present investigation was to assist NASA in 
theoretically evaluating the design of the model pilot quiet tunnel through 
implementation of the UARL boundary layer prediction procedure. The first step 
in fulfilling this purpose was to evaluate the UARL procedure's ability to 
predict the behavior of boundary layers subjected to large favorable pressure 
gradients. Although the ability of the UARL procedure to predict the boundary 
layer mean flow behavior in going from the laminar to the turbulent state had 
been previously extensively assessed (refs. l4 and 15 )> no detailed assessment 
had been previously made of the procedure's ability to predict the effects of 
large favorable pressure gradients on an initially turbulent boundary layer. 

The experimental highly accelerated boundary layer data available for compari- 
son with theoretical predictions falls into three broad areas ; adiabatic wall 
low-speed flows, adiabatic wall moderately supersonic flows, and low-speed flows 
in the presence of wall heat transfer . The amount of data available with heat 
transfer in highly accelerated turbulent flows, in sufficient detail to enable 
predictions to be made, is limited. However, enough data is available to make 
some meaningful comparisons. 


The Effect of Pressure Gradient Upon Transition 

In general, it is expected that relaminarization will occur in regions of 
strong favorable pressure gradient and as a first step in evaluating the ability 
of the model to predict relaminarization, a series of calculations were made to 
compare theoretical predictions of the forward transition Reynolds number with 
the experimental data of Feindt (ref. 17) over a range of pressure gradients. 

The results of several calculations made for various pressure gradients are 
shown in Fig. 1 where the predicted variation of the momentum thickness Reynolds 
number at transition, Re^ with the pressure gradient parameter (0‘^/v)dU/dx, is 
compared with the theory of van Driest and Blumer (ref. l 8 ) and with the 
experimental data of Feindt (ref. 17) as presented by Hall and Gibbings (ref. 
19 ). The predicted transition location was taken to be at the point of minimum 
skin friction. The UARL predictions and the data of Feindt are for a free- 
stream turbulence level of 1.2 percent. The van Driest -Blumer theory is 
presented for 1 percent free -stream turbulence. Both theories, as well as the 
experimental data, show only a mild effect of pressure gradient on transition 
Reynolds number over this moderately low range of the pressure gradient param- 
eter. The UARL predictions are in good agreement with Feindt 's data, showing 
an increase in transition Reynolds number with increasing pressure gradient 
parameter . 


10 



The Effect of Large Strearawise Accelerations 
Upon the Boundary Layer in Supersonic Nozzles 

As part of the present investigation predictions of the effect of large 
streamwise accelerations upon initially fully-turbulent boundary layers were 
made corresponding to boundary layer measurements obtained experimentally by 
Nash-Webber (ref. 6). The measured boundary layers were developed in nozzles 
constructed with a flat plate as one wall and a contoured surface as the 
opposite wall; the experimental data were taken for the boundary layer develop- 
ing on the flat plate wall. The nozzles, referred to as nozzle "A" and nozzle 
"c" in ref. 6, had exit Mach numbers of approximately 2.0 and had adiabatic 
walls. Both nozzles subjected the boundary layers to favorable pressure 
gradients. Since Nash-Webber (ref. 6) did not report a free-stream 
turbulence intensity, a level of 2 percent was assumed at the nozzle inlet. 

The free-stream turbulence kinetic energy, is assumed to remain constant 
under either acceleration or deceleration, therefore, the intensity 
(^/3)t/Ue, varies inversely with free-stream velocity. Comparisons between 
predictions made using the UARL transition model, fully-turbulent theory 
(frozen mixing length), and the data of Nash-Webber (ref. 6) for three different 
stagnation pressures in nozzle "A" are shown in Figs. 2 through 4. Variation 
of the stagnation pressure provides a means of varying the free-stream unit 
Reynolds number, Reynolds number increasing with increasing pressure, while 
maintaining the same pressure gradient. The acceleration parameter, K, 

K = rfP/tf. (1) 


however, is reduced as the stagnation pressure is increased due to the increase 
in density. The results shown in Figs. 2 through 4 indicate that the predicted 
variation of momentum thickness Reynolds number and shape factor made using the 
UARL prediction procedure are in good agreement with the experimental data for 
the three values of stagnation pressure considered. In Figs. 2 through 4 (as 
well as Figs. 5 and 6) the distance along the nozzle has been nondimensionalized 
by a reference length, Lygf = 30.5 cm (12 in.). The variation of the dimen- 
sionless mixing length (jtoo/6), indicates that the boundary layers effectively 
relaminarize at a dimensionless distance of about 2.0. The r elaminarization 
is also indicated by the sharp increase in the boundary layer shape factor. 

The predictions made using fully-turbulent theory are significantly inferior 
as they do not correlate well either with the experimental data or the 
predictions made using the transition model past the region of relaminarization . 
The predictions ofthe boundary layer development in nozzle "C" at two values 
of stagnation pressure are shown in Figs. 5 and 6. The predictions made using 
the UARL procedure are again in good agreement with the experimental data for 
both values of stagnation pressure, however, the strearawise variation of the 
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momentum thickness Reynolds number at a stagnation pressure of 1.66 x 10"! atm, 
shown in Fig. 5, is somewhat overpredicted. The sharp drop in the dimensionless 
mixing length occurs at approximately x= 2.75 for both values of stagnation 
pressure indicating the region of relaminarization, and it is downstream of 
this region that the fully-turbulent predictions deviate significantly frcrn 
both the data and the transitional predictions. 


The Effect of Large Streamwise Accelerations 
Upon Low -Speed Boundary Layers 

Rredictions for strongly accelerated flows. - A detailed investigation of 
a low -speed highly-accelerated boundary layer which is suitable for ccmparison 
with theory was carried out by Blackwelder and Kovasznay (ref. 20). In this 
investigation a flow in which the free-stream velocity was initially 295 cm per 
sec was accelerated smoothly to a velocity of 1260 cm per sec over a distance 
of 350 cm. During -toe acceleration the acceleration parameter, K, reached a 
maximum of U.8 x 10”® which is well above the usual value at which relaminari- 
zation is observed, niis naximum acceleration parameter is reached at 93^ cm 
and the acceleration then decreases to zero at 1200 cm. The free-stream 
turbulence level was of the order of 0.5 percent as determined from measurements 
of u', V*, and w'. Theoretical predictions of the flow were made both with the 
UARL transition procedure using 0.5 percent free-stream turbulence and with the 
fully-turbulent (constant dimensionless mixing length) theory. Comparison of 
the two predictions with the experimental data axe presented in Figs. 7 and 8. 
Figure 7 shews the variation of the momentum thickness Reynolds number and the 
shape factor. The predictions made using the UARL procedure are in good agree- 
ment with the experimental data. In contrast with the transition prediction, 
the fully-turbulent prediction does not give good quantitative predictions of 
the manentum thickness Reynolds number and shape factor downstream of the 
point of maximum acceleration. The variation of the skin friction coefficient 
and the dimensionless mixing length is shown in Fig. 8. The prediction of skin 
friction coefficient made using the UARL procedure shows qualitative agreement 
with the data but quantitative disagreement. Both theory and experiment show 
a drop in skin friction followed by a rise, however, the theoretically 
predicted rise anticipates the experimentally observed rise. It should be noted 
that the fully-turbulent theory is not even in qualitative agreement with 
experiment as far as skin friction or shape factor predictions are concerned. 

Tbe dimensionless mixing length shows that relaminarization is effectively 
complete by x = 9^0 cm. This is close to the point of maximum acceleration 
and coincides with the point where the fully-turbulent prediction begins to 
show large discrepancies with the data. The dimensionless mixing length then 
shows a sharp increase Indicating that transition back to turbulent flow occurs 
in the downstream section where the flow is subjected to zero pressure gradient. 
This result is consistent with the experimental results of Blackwelder and 
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Kovasznay (ref. 20). With regard to the predicted and experimental values of 
the skin friction coefficient it is of interest to note that the boundary- 
layer is initially a f ully-turbulent , zero pressure gradient boundary layer. 

At the initial value of the momentum thickness Reynolds number of 2500, the 
data correlation of Smith and Walker (ref. 21 ) would indicate that the skin 
friction coefficient should be close to 0.00315 > a value in close agreement 
with the predicted value, but higher than that measured by Blackwelder and 
Kovasznay. Thus, some error in the experimental values of skin friction, which 
were determined by measurement of the velocity gradient in the viscous sublayer 
near the wall, may be present. 


Predictions for sink -flows. - Jones and Launder (ref. 22) have made 
experimental measurements of sink-flow turbulent boundary layers which are of 
particular interest because they are the only turbulent boundaiy layers with 
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length scales may develop at the same rate (ref. 23). In sink -flows the 
acceleration parameter, K, is constant and the momentum thickness Reynolds 
number and shape factor obtain approximately constant asyn^itotic values over a 
considerable strearawise distance. The velocity profiles also become similar 
once these asyn^totic conditions are reached. Two predictions were made under 
conditions investigated by Jones and Launder (ref. 22) for K = 1.5 x 10"^ and 
K = 2.5 X 10"^. The theoretical predictions for a given value of K were carried 
out until a sensibly constant value of the momentum thickness Reynolds number 
was reached. Curves and experimental data taken from ref. 22 for the asymptotic 
values of momentum thickness Reynolds number and shape factor as a function of 
the acceleration parameter are shown in Fig. 9. Th® symbols represent data 
taken from five sources (refs. 24 through 28), all of which show relaminariza- 
tion except for tbe data of Julien, et al. (ref. 26) taken in the Stanford 
tunnel. This discrepancy has been noted by Jones and Launder (ref. 22) who 
attribute it to the inability of the Stanford tunnel to produce high-K sink- 
flows because of no independent means of adjusting Resmolds number. However, 
an alternative possibility may be the presence of some turbulence generating 
mechanism. Because of the discrepancy, no comparisons between theoretical 
predictions and the Stanford data for high accelerating flows were made except 
for the Moretti heat transfer runs discussed subsequently. The triangular data 
points indicate the new data of Jones and Launder (ref. 22) and the flagged 
circles represent the asymptotic values obtained from the UARL prediction 
procedure. The theoretically predicted and the experimentally measured values 
of Rg and H are in excellent agreement. A more stringent test of the theory is 
the comparison of the predicted and measured velocity profiles when the 
momentum thickness Reynolds number becomes constant. The predicted velocity 
profile for K = 1.5 x 10“^ is compared with the experimental data at three 
values of R0 for the given value of K in Fig. 10. The agreement between the 
prediction and the data is quite good, particularly for the experimental pro- 
file at R9 = 711. Some discrepancy is evident for the other -two experi- 
mental profiles in the vicinity of a normal distance Reynolds number of one. 
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Similar results are presented in Fig. 11 for K = 2.5 x 10"^. The comparisons 
are generally good at this higher value of K with discrepancies occurring again 
at a normal distance Reynolds number of one. The comparisons shown in Figs. 

10 and 11 indicate that the UARL prediction procedure may not adequately allow 
for the thickening of the viscous sublayer, which is known to occur during 
strong acceleration. 

Effect of acceleration on heat transfer. - Three comparisons between 
theoretical predictions and experimental data were made to determine how well 
the UARL procedure could predict the variation of Stanton number in a strongly 
accelerated flow. The three experimental cases which were used for comparison 
were measured by Moretti and Kays (ref. 29) and termed Runs 12, 4l, and F-2. 

All three of these experimental cases were basically similar in that a constant 
velocity low -speed flow was first subject to a nearly step reduction in wall 
ten^jerature . The flow was then allowed to undergo some thermal development 
before accelerations of varying strength were applied. In Run 12, the 
acceleration parameter increases sharply to a maximum K of 3*39 x 10"^ and 
then rapidly returns to zero. The comparison between the predicted and measured 
variation of Stanton number along with the wall temperature variation for Run 12 
is shown in Fig. 12. The location of the initiation and end of the acceleration 
are indicated in the figure. Both the prediction using the UARL transition 
model and the experimental data indicate that there is no effect on the Stanton 
number until the dimensionless distance x = 5*0 is reached with good agreement 
to that point even though the acceleration was initiated at x = X = 

5.0 the UARL prediction and the data indicate a sharp decline in Stanton number 
although the prediction gives a larger decrease than that experimentally 
observed. The fully-turbulent prediction shows no decrease. At x = 5*5 the 
data show an increase in Stanton number. This behavior also is predicted by 
the transitional theory, however, the predicted increase is somewhat delayed 
compared to the data. The variation of the dimensionless mixing length, also 
shown in Fig. 12, shows that effectively relaminarization occurs just past 
X = 5«0 and also indicates that the failure of the UARL prediction to follow 
the data once the pressure gradient is removed may be due to incomplete 
transition back to turbulent flow. However, on a qualitative basis, the 
prediction made using the UARL procedure shows significantly better results 
than the prediction made using fully-turbulent theory. Figure 13 shows similar 
results for Run kl. In Run 4l the maximum value of K is 1.8U x 10 and K is 
nearly constant at this value over the distance from x= ^*0 to x= 5«0. 

Again, qualitative agreement is shown between the data and the UARL predictions. 
However, although the data only indicate a slight degree of relaminarization, 
the theory indicates a strong relaminarization. Based upon the sink -flow 
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boundary layer results of Pig. 9, a flow subject to an acceleration parameter 
of 1.8 X 10~^ should show a significant relaminarization effect. Figure lU 
then shows the results for Run F-2 where K reached a maximum value of 3 *0^ x 
10~° and was nearly constant from x= ^*0 to^ X= 5»3* Here, too, the 
prediction using the transition model shows qualitative agreement. From the 
results shown in Figs. 12 through l4 it is noted that although only qualitative 
predictions of Stanton number through the strongly accelerated regions are 
obtainable, the onset of the effect of acceleration is well predicted by the 
UARL procedure, whereas the fully-turbulent theory gives no indication of 
feeling the effect of acceleration. 
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CALCULATIONS FOR A QUIET TUNNEL 


The quiet tunnel program currently being undertaken by NASA (ref , -l6) mses 
a converging -diverging nozzle with an annular slot just upstream of the nozzle 
throat. The slot is used to remove the boundary layer which develops on the ' 
inlet duct and convergent nozzle section walls thus reducing the boundary 
layer thickness at the throat station. Predictions of the boundary layer 
behavior in the inlet duct and approach section of the NASA nozzle (convergent 
portion) were made to give an indication of the nozzle slot height required to 
remove the approaching boundary layer. A new boundary layer then develops 
from the edge of the slot and continues to grow as the flow moves downstream 
into the supersonic section of the nozzle. If the tunnel wall boundary layer 
in the supersonic portion of the nozzle undergoes transition and becomes 
turbulent, the acoustic radiation emanating from the boundary layer will 
dominate the process occurring oii a model plac^ in the tinpel and ^ 

result in the measurements being tunnel dependent. !fhj^efore, 

it is nece^sse^ to aaintain a laminar boundary layer in the supersonic section 
of the nozzle as far downstreeun as possible. With these considerations in 
mind, a parametric study of the effects of free-stream turbulence, heat 
transfer, Reynolds number, acceleration and Mach number on the boundary layer 
development in the supersonic section was conducted to determine both if the 
boundary layer developing in supersonic regions of the contemplated NASA tunnel 
would remain laminar and how varying various parameters would effect the 
transition location. As a portion of these parametric studies, inviscid flow 
pressure distributions were computed for various nozzle geometries to facilitate 
the calculations of Mach number and acceleration effects. In addition, 
transonic and supersonic inviscid flow computations were made for the nozzle 
geometry specified by NASA to con 5 >are with the pressure distribution supplied 
by NASA in these regions of the nozzle. 


Inviscid Flow Calculations 

Pressure distribution in the transonic section of the NASA nozzle. - 
Although pressure distributions were supplied by NASA for the transonic and 
supersonic sections of the NASA nozzle, several inviscid flow field predictions 
were made for NASA specified nozzle geometry using United Aircraft Corporation 
inviscid computer codes. The calculations were made at NASA's request to 
compare the inviscid flow predictions made using UAC codes with the predictions 
made using NASA codes. The nozzle coordinates and pressure distribution 
supplied by NASA are given In TABLE I . Transonic flow field conqoutations were 
performed using the UAC SUB code (ref. 30). The present version of the SUB 
code cannot account directly for the presence of the slot in the nozzle wall, 
thus, two predictions of the transonic flow region were made each modifying 


16 





the approach section nozzle wall to represent the slot. Since the computation 
of the inviscid flow field in the supersonic section of the nozzle is. dependent 
on the Msich number and pressure distributions predicted from the SUB code in 
the region of the sonic line, the location of the resulting sonic lines from 
the two predictions were then compared to determine how sensitive the sonic line 
was to changes in the upstream geometry. 

Ihe first prediction was made for a nozzle in which the slot lip was 
smoothly fared into the approach section nozzle wall. Ihe predicted stream- 
lines and sonic line for this contour, referred to as contour "A”, are shewn 
in Fig. 1^. Ihe broken line indicates the original NASA contour and the heavy 
solid line indicates the c<mitour used in the computation. Ihe NASA contour and 
contour "A" are identical downstream of the slot, llie second computation was 
made for an approach section outer contour which corresponds to the estimated 
steignation strecuoline which divides the flow passing through the nozzle throat 
from that pasat^ throuii^ the slot. l!he coordinates of this streamline were 
obtained by mass <^servation of the flow throu^ the sonic slot. Hae predicted 
streamlines and sbeSc lii» for this contour, referred to as contour "b", are 
shown in Fig. l6. Downstream of the slot contour "B" is also Identical to the 
NASA contour. An enlarged cooq>arison of the sonic lines predicted using 
contours "A" and "B” presented in Fig. 17 shows that the effect of the modifi- 
cations to the nozzle wall does not affect significantly either the shape or 
location of the sonic line. A ccmqMrison between the wall pressure distributi(»i, 
normalized by stagnation pressure, predicted using contours "A" and "b" and tne 
NASA distribution is shown in Fig. I8. ibe pressure distributions for contours 
"a" and ”B", and the NASA pressure distributicxi show only small differences in 
the upstream region between x= "3«0 and x= and the distributirais are 
nearly Identical in the throat region between x= “0*5 and x = 0.0. The 
pressure distributicxis downstream of the throat predicted for both contours 
"A" and "B" are slightly different from the NASA distribution. Furthermore, 
the NASA distribution shows a small discontinuity in slope at x = 0.350 which 
does not appear in either UAC calculation. It should be pointed out that the 
NASA pressure distribution was calculated using the theory of Hopkins and Hill 
(ref. 31) from x= “O.90 to x= O.29U and using the method of characteristics 
downstream of Xf 0,294. Ibus, the apparent discontinuity of slope in the NASA 
pressure dist^lbc^ion may be the result of slightly incompatible Initial 
conditions for the method of cheuracteristics computatim. 

Pressure distributions in the sunersoiic section of the nozzle. - A 
computation of the flow field in the stq>ersonlc section of the NASA nozzle was 
made using the UAC SUIER code (ref. 32) with initial Mach number and pressure 
distributions obtained from the SUB code along a plane pe3rpendicular to the 
nozzle axis at x= 0.25. n»e predicted wall pressure distribution, shown in 
Fig. 19» is in excellent agreement with the NASA pressure distribution. Ibe 
predicted flow field contains no e3q>ansion fans or shocks and the predicted 
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exit plane Mach number is close to uniform with a maximum of U.98 at the 
centerline and a minimum of 4.93 at the nozzle wall. Figure 20 shows an 
enlarged comparison between pressure distributions in the vicinity of the 
throat as predicted by the SUB and SUIER codes, and that specified by NASA. 

All three pressure distributions show good agreement. Ihe SUIER prediction 
fares smoothly into the prediction made using the SUB code at x = 0.25 . The 
SUEER prediction Initially gives slightly higher values of wall pressure ratio 
than the NASA prediction, however, by X= 0.5 both distributions are in good 
agreement. Neither the SUB prediction, nor the combined SUB-SUIER prediction 
exhibit the discontinuities of slope observed in the NASA pressure distribution. 

Invisold calculations for acceleration and Mach number effects. - Several - 
additional Inviscid flow field calculations were made using the UAC Perfect 
Nozzle code to obtain pressure distributions to be used in boundary layer 
calculations showing the effect of acceleration and the effect of Mach number. 
The Perfect Nozzle code generates the coordinates and pressure distribution of 
a sharp corner throat bell nozzle with a uniform exit Ifech number for a 
specified throat -to -exit area ratio. Nozzle contours and pressure distributions 
were computed for area ratios corresponding to exit Mach numbers of 4.1, 5.0, 
and 6.0. The resulting wall pressure distributions for these nozzles are 
shown in Fig . 21 . 


Viscous Flow Calculations 

Boundary layers in subsonic portion of NASA nozzle. - Calculations of the 
boundary layer flow field on the inlet duct and nozzle approach section walls 
of the NASA nozzle were carried out using the pressure distribution specified 
by NASA given in TABLE II, with a wall temperature of 65.5 deg C (150 deg F), 
a free-stream ten5>erature of 104.4 deg C (220 deg F), and a total pressure of 
10.2 atm (150 psia). The inlet duct extends from x= "6.91 to x= -0.9II and 
the approach section extends from x= “O.911 to x = -0.023 where x = */^ef 
I^ef is taken as 30.5 cm (12 in.). Calculations performed for free-stream 
turbulence levels of 0.I8 and 1.0 percent based on a reference velocity of 
approximately II76 cm per sec ( 4 o ft per sec) were initiated from an assumed 
stagnation point at the upstream end of the inlet duct. These turbulence 
levels correspond to 1 and 6 percent based on the actual initial velocity in 
the inlet duct of 213 cm per sec (7 ft per sec). Transition from laminar to 
turbulent flow was predicted to occur in the inlet duct well upstream of the 
approach section of the nozzle for both turbulence levels. The predicted 
boundary layer behavior in the approach section of the nozzle downstream of the 
inlet duct is shown in Figs. 22 through 24 where the boundary layer thicknesses 
are nondimensionalized with respect to the nozzle throat radius of 1.27 cm 
(0.5 in.). In the nozzle approach section the boundary layer thickness 
decreases continuously under the influence of the strong favorable pressure 
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gradient for both turbulence levels, as shown in Fig. 22. The distribution of 
the dimensionless mixing length, also shown in Fig. 22, indicates that for 
1.0 percent free-stream turbulence relaminarization of the boundary layer is 
initiated at x = -0.65 and is nearly complete at x= The mixing length 

distribution predicted for 6.0 percent free-stream turbulence indicates that 
relaminarization is completed over a longer distance being initiated at x = 

-0.6? and showing effectively complete relaminarization by x= -0.3. However, 
for both turbulence levels the boundary layer is laminar at , the throat. The 
momentum thickness also decreases continuously for both levels of turbulence 
as the throat is approached, as shown in Fig. 23. The behavior of the boundary 
layer thickness near the throat of the nozzle is shown on an enlarged scale in 
Fig. 2U with an indication of the location and height of the slot. As can be 
seen in this figure, the predicted boundary layer thickness is almost identical 
to the slot height, thus indicating that the slot is large enough to remove the 
entire velocity boundary layer. 

Boundary layer downstream of the slot in the NASA nozzle. - Predictions 
of the boundary layer development in the transonic and supersonic section of 
the NASA nozzle were made for free-stream turbulence levels of 1, 3, and 6 
percent based on a reference velocity of 1176 cm per sec (40 ft per sec). The 
results of these calculations, which were initiated at the slot with a 
negligible displacement thickness (6* = 0.000254 cm = 0.0001 in.), are presented 
in Figs. 25 through 27. The flow conditions are indicated on the figures. In 
Figs. 25 through 27, as well as all the figures which follow, the boundary 
layer thicknesses are normalized by the nozzle throat radius r-^ = 1.27 cm 
(0.5 in.) and axial distances are nondimensionalized by the reference length 
Lref = 30.5 cm (12 in.). The acceleration parameter, shown in Fig. 25, has a 
maximum value of 1.53 x 10“^ at the slot and decreases continuously through the 
throat. At an axial distance of x= 0.1, the acceleration parameter has 
dropped to less than 10 percent of its maximum value and between x= 0.1 and 
the exit of the nozzle, at x= 1*289, K decreases almost linearly to a value 

of 0.03 X 10 “^. 

The incon 5 )ressible shape factor, H^, i.e., the locally evaluated shape 
factor obtained from the incoinpresslble definitions of momentum and displace- 
ment thicknesses, is also presented in Fig. 25. The incompressible shape 
factor is presented rather than the true shape factor since it more satisfac- 
torily indicates the laminar or turbulent state of the boundary layer when 
density gradients across the boundary layer are significant. For the 1 percent 
turbulence level H decreases sharply under the influence of the strong 
acceleration in the throat region. As the acceleration is relaxed, the 
incompressible shape factor increases to a relatively constant value of about 
2.5 indicating a laminar boundary layer for the entire nozzle length. For the 
6 percent turbulence level H* shows the same initial trend, but near x= 0*1 
H* drops to a constant value of 1.4 which indicates that the boundary layer has 


19 


undergone transition to turbulent flow. This behavior is also observed at 
X = 0.2 for 3 percent free -stream turbulence. The streainwise variations of the 
momentum thickness and the compressible shape factor are shown for the three 
levels of free -stream turbulence in Fig. 26. For the higher free-stream 
turbulence levels the turbulent nature of the boundary layer is indicated by 
the more rapid grcwth of the momentum thickness. Finally, the variations of 
the dimensionless mixing and boundary layer thicknesses are shown in Fig. 27« 
The variation of the dimensionless mixing length shows that the boundary layer 
remains laminar for the 1 percent turbulence level. The sharp rise to nearly 
constant values of 0.11 for the higher turbulence levels Indicates the rapid 
transition to a turbulent boundary layer in the presence of the higher free- 
stream turbulence . 


Parametric Variation in the Supersonic Portion of the Nozzle 

Effect of heat transfer. - Since wall cooling is known to have a stabili- 
zing effect for a constant pressure boundary layer, it may also be expected to 
maintain a laminar boundary layer in a nozzle where large favorable 
pressure gradients are present. Thus, predictions of the effect of wall 
cooling on the boundary layer in the NASA nozzle were made for wall temperatures 
of -17.8 deg C (O deg F) and -73 «3 deg C (-100 deg F) for turbulence levels of 
1 and 3 percent. The results of the calculations made for 1 percent free- 
stream turbulence are compared with predictions made for an adiabatic wall in 
Figs. 28 through 30. These cooqparisons indicate that for this low value of 
free-stream turbulence wall cooling promotes earlier boundary layer transit ic»i 
as can be seen in the plots of incompressible shape factor, presented in Fig. 

28, momentum thickness, shown in Fig. 29, and dimensionless mixing length, 
shown in Fig. 30. Predictions made for the same wall tenperatures and flow 
conditions with a free-stream turbulence level of 3 percent are presented in 
Figs. 31 through 33. The predicted variations of incoopressible shape factor, 
momentum thickness, and mixing length indicate that for this higher turbulence 
level wall cooling has little effect on the location of transition. 

The apparent discrepeincy of the prediction of wall cooling promoting 
transition can be explained by citing the dual role of wall cooling in the 
presence of a pressure gradient. In the absence of wall roughness and pressure 
gradients, wall cooling stabilizes the boundary layer retarding transition. 
However, in the presence of a favorable pressure gradient this direct stabiliz- 
ing effect is offset by the destabilizing effect resulting frcan the reduction 
of the acceleration parameter, K, as is shown in Figs. 28 and 31. This reduc- 
tion of the acceleration paraimeter is due solely to the reduction of the 
kinematic viscosity near the wall. It appears that this dual effect of wall 
cooling may either enhance or inhibit transition dej)ending upon whether the 
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direct stabilizing effect of wall cooling is greater or less than the indirect 
destabilizing effect of the reduced acceleration parameter. Since the 
transition point varied with wall temperature for the case of 1 percent free- 
stream turbulence while showing relative insensitivity to wall temperature at 
a 3 percent free -stream turbulence, the integrated effects of wall cooling 
could be, and apparently were, different for the two different turbulence 
levels. The important observation is that in the presence of a favorable 
pressure gradient wall cooling may not necessarily inhibit transition. 

Effect of Reynolds number, - The effect of variation of the free -stream 
unit Reynolds number on the boundary layer development in the transonic and 
supersonic sections of the NASA nozzle was achieved by varying the nozzle 
stagnation pressure; the variation of stagnation pressure then being reflected 
in the value of the kinematic viscosity. The results of the calculations, which 
were performed for stagnation pressures of 5-1 j 10.2, and 25.5 atm (75j 150, and 
375 psia), are shown in Figs. 3^ through 36. All calculations were made for a 
free-stream turbulence level of 1 percent. The calculation made using a 
stagnation pressure of 5*1 atm (75 psia) corresponds to a low Reynolds number 
case, whereas the 25.5 atm (375 psia) stagnation pressure calculation represents 
a higher Reynolds number. Within this range of stagnation pressure the free- 
stream unit Reynolds number at a given axial location varies by a factor of 5 
between the three calculations. The variation of the incompressible shape 
factor, shown in Fig. 3^> indicates that for both the low ( 5 .I atm (75 psia)) 
and intermediate (10.2 atm (150 psia)) values of Reynolds number a laminar 
boundary layer is maintained for the entire length of the nozzle; however, the 
higher Reynolds number prediction shows that the boundary layer undergoes 
transition from laminar to turbulent slightly downstream of the throat. The 
acceleration parameter, also shown in Fig. 3^» shows substantial variation with 
Reynolds number. This is a direct result of the variation of the kinematic 
viscosity. The variation of the momentum thickness and compressible shape 
factor along the nozzle are shown in Fig. 35* These results again indicate 
the laminar state of the boundary layer for both low and intermediate Reynolds 
numbers. It is also noted that the predicted variation of the momentum thick- 
ness with axial distance for low and intermediate Reynolds numbers, and for 
the higher Reynolds number prediction prior to transition, is consistent with 
the usual laminar boundary layer theory; the momentum thickness grows more 
rapidly for lower Reynolds numbers. Finally, the variation of the dimension- 
less mixing length and the boundary layer thickness is shown in Fig. 36 . The 
variation of the dimensionless mixing length indicates that transition occurs 
only for the higher Reynolds number prediction. 

Effect of acceleration. - Calculations were made to determine the effect 
of varying the acceleration parameter upon the boundary layer development. The 
calculations were made using the pressure distribution obtained from the 
Perfect Nozzle deck for a Mach 5 exit sharp corner throat, shown in Fig. 21. 

The sharp corner throat represents a zero radius of curvature and gives the 
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greatest possible wall pressure gradient in the throat region. The resulting 
nozzle geometry generated by the UAC Perfect Nozzle deck differed from the NASA 
nozzle primarily in the region of the throat leading to different static 
pressure distributions. The effects of this variation of the nozzle geometry 
on the acceleration parameter, incompressible shape factor, and momentum 
thickness are shown in Fig. 37 where they are compared with the values 
predicted for the specified NASA nozzle geometry. The calculations for the 
sharp corner throat nozzle were initiated at the throat with a displacement 
thickness equal to that at the throat of the NASA nozzle. The acceleration 
parameter for the sharp cornered throat nozzle only differed slightly from 
that of the NASA nozzle, as shown in Fig. 37 • The slight difference in 
acceleration parameter did not lead to any significant difference in boundary 
layer development, as shown in Figs. 37 and 38 * The major differences between 
the two calculations were in the throat region where differences were predicted 
in the incompressible shape factor. 

Fffect of exit Mach number. - The effect of variation of the nozzle exit 
Mach number was studied by making boundary layer predictions for the pressure 
distributions shown in Fig. 21 for exit Mach numbers of 4.1, 5, and 6. The 
pressure distributions are for nozzles having sharp corner throats and 
uniform exit flow. The predictions were initiated at the nozzle throat with a 
displacement thickness equal to that predicted for the NASA nozzle with a 
free-stream turbulence of 1 percent. The predicted variations of acceleration 
parameter, incompressible shape factor, and momentum thickness are compared in 
Fig. 39 * The variation of the acceleration parameter shows that there is very 
little effect due to exit Mach number in the throat region; however, the 
characteristic linear decrease in acceleration parameter downstream of the 
throat shows that as the exit Mach number is increased, the acceleration 
parameter decreases more slowly. The variation of the incompressible shape 
factor shows a drop near the throat and then an increase to a relatively 
constant value of 2.5 in the Mach 5*0 and Mach 6.0 nozzles while continually 
increasing to a value of 2.7 in the Nfech 4.1 nozzle. These variations of the 
incompressible shape factor, as well as the smooth continuous rise in momentum 
thickness, indicate that a laminar boundary layer is maintained for the 
entire length of all three nozzles. The variations of the compressible shape 
factor, dimensionless mixing length, and boundary layer thickness are shown 
in Fig. 40. The variation of the dimensionless mixing length indicates the 
laminar state of the predicted boundary layers. It is interesting to note that, 
although still small, the value of the dimensionless mixing length at the nozzle 
exit shows a trend of increasing with exit Mach number. This trend is probably 
due, at least in part, to the longer distance over which the boundary layer 
develops in higher Mach number nozzles. 
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DISCUSSION 


The inconsistencies found when attempting to correlate transition data from 
several supersonic wind tunnels (refs. 3 and 4) and the evidence that the cause 
of these inconsistencies is the aerodynamic noise radiating from the tunnel wall 
turbulent boundary layer clearly indicates the need for a quiet supersonic wind 
tunnel. In addition, a recent investigation by McDonald and Kreskovsky (ref. 

33) shows that the level of free -stream turbulence can significantly alter the 
structure of even fully-turbulent boundary layers and since the fluctuating 
pressure associated with an acoustic disturbance can be easily related to a 
fluctuating velocity (ref. 15), it is reasonable to assume that the varying 
magnitude of aerodynamic noise found in different tunnels could effect the 
fully-turbulent boundary layer development as well as transition location. 

Tile re Tore , to aid in the design of a quiet tunnel, a parametric study of boundary 
layers in supersonic nozzles has been conducted in this investigation. 


Prior to performing the parametric study it was necessary to ascertain 
that the boundary layer prediction procedure used would give accurate predictions 
of relaminarizing flows. Comparisons of theoretical predictions of momentum 
thickness and shape factor with the supersonic data of Nash-Webber (ref. 6) and 
the low-speed experimental data of Blackwelder and Kovasznay (ref. 20) were 
very good. However, comparison between theoretical predictions of skin friction 
and the data of ref. 20, although showing qualitative agreement, showed 
quantitative discrepancies. Similarly, comparisons between theory and experi- 
ment for sink-flow boundary layer shape factors and momentum thicknesses were 
very good, however, comparisons of velocity profiles, although very good in the 
outer region, showed discrepancies in the wall region. Finally, predictions 
of Stanton number through relaminarization (when compared with the data of 
Moretti and Kays (ref. 29)) indicate that the onset of the effect of acceleration 
is accurately predicted although the actual variation of Stanton number through 
relaminarization is somewhat overpredicted. Thus, in general, the theory 
predicts the initiation of relaminarization and the variation of integral 
thicknesses through relaminarization very well. Comparisons between theory and 
experiment for skin friction and Stanton number, although showing qualitative 
agreement, show quantitative discrepancies. However, as discussed in the 
Comparison section, neither skin friction nor Stanton number data are 
completely reliable and in any event present a much more difficult measurement 
task than the mean flow or integral thicknesses. 

The subsequent calculations made for the quiet tunnel study indicate that, 
perhaps not surprisingly, transition from laminar to turbulent flow is sensi- 
tive to the level of free-stream disturbance (see Figs. 25 through 27); a 
laminar boundary layer being maintained only for the lowest turbulence level 
prediction. With regard to the effect of wall cooling, it was found (see Figs. 
28 through 33) that a destabilizing effect may result from wall cooling in the 
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presence of favorable pressure gradients. A possible explanation for this 
result is the dual role played by -wall cooling in the presence of a favorable 
pressure gradient, the direct stabilizing effect and the destabilizing effect 
of reduction of the acceleration parameter. Predictions of the effect of 
Reynolds number indicate that, as one would expect, increasing Reynolds number 
promotes transition (Figs. 3^ through 36). Predictions of the effect of 
variation of the pressure gradient (via variation of nozzle geometry) show 
only slight differences for the geometries examined (see Figs. 37 and 38). A 
significant variation of the pressure distribution in the supersonic section of 
the nozzle cannot be achieved since the major portion of the potential energy 
associated with the flow is converted to kinetic energy in the transonic 
region near the throat. Finally, calculations made for the effect of Mach 
number show that laminar boundary layers may be maintained for the entire 
nozzle length for nozzles with exit Mach number of at least 6.0 for low free- 
stream turbulence levels (see Figs. 39 nnd 4o). 


CONCLUSIONS 


Conclusions Regarding the Assessment of the 
UARL Boundary Layer Prediction Procedure 

1. Hie UARL boundary layer prediction procedure, previously shown by 
McDonald and Fish to account for the influence of such items as free-stream 
turbulence and wall roughness upon the transition from laminar to turbulent 
flow in a zero pressure gradient boundary layer, has been shown in the present 
study also to account for the movement of the transition region due to 
moderate pressure gradients. 


2. The UARL procedure was also found to accurately predict the 
variation of boundary layer integral paranssters in an initially turbulent 
boundary layer subject to large accelerations in both low -speed and supersonic 
flows . 


3. Hie UARL procedure apparently predicts the onset of the effect of 
strong acceleration on wall -dominated properties such as Stanton number. 

However, the variation of wall -dominated properties in boundary layers under- 
going relaminarization is in only qualitative agreement with presently available 
data. 


Conclusions Regarding the Quiet Tunnel 

1. Ihe predicted boundary layer in the approach section of the NASA 
nozzle effectively relaminarizes for free-stream turbulence levels as high as 
6 percent based upon an inlet velocity of 213 cm per sec (7 ft per sec). 

2. The predicted boundary layer thickness at the slot is approximately 
equal to the slot height. 

3. It should be possible to maintain a laminar boundary layer for the 
full length of the NASA nozzle with adiabatic walls, provided that the free- 
stream turbulence level is of the order of 6 percent or less in the approach 
section based upon a mean velocity of 213 cm per sec (7 ft per sec). 

k. The predicted effect of wall cooling promotes early transition at 
1 percent free-stream turbulence but shows little effect on transition 
location at 3 percent free-stream turbulence. 
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5. Transition is predicted to occur earlier with higher free -stream 
unit Resmolds number. 

6. Variation of the acceleration parameter downstream of the nozzle 
throat could not be increased significantly by changes in geometry and only 
insignificant effects of the resulting changes in acceleration were observed 
on the boundary layer . 

7. The indications are that a laminar boundary layer may be maintained 
in nozzles with exit Mach numbers as high as 6.0. 


26 


APPENDIX A 


The Basic Equations 

Within the framework of boundary layer theory, various authors (for 
example, Schubauer and Tchen (ref. 3^)) have reduced the time-averaged Navier- 
Stokes equations to the compressible boundary layer equations of motion. For 
two-dimensional or axisymmetric flows, steady in the mean, the boundary layer 
approximations to the momentum, energy, and continuity equation become 
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where x and y are coordinates in the streamwise and transverse directions, 
u and V are velocity components in the x and y directions, p is density, P is 
pressure, Cp is specific heat, T° is total temperature, r is the radius of 
curvature for an axisymmetric body, and the exponent a is zero for two-dimen- 
sional flows and unity for axisymmetric flows. The effective shear stress, t, 
and the effective heat transfer, Q, are defined as 
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where u is viscosity, k is thermal conductivity, T is static temperature. In 
Eqs. (A-l) through (A - 5 ) overbars indicate averaged quantities and primes 
indicate fluctuating quantities. The equations are valid for laminar, 
transitional, or turbulent flows; obviously, for laminar flows the primed 
quantities are zero. For turbulent and transitional boundary layers it is 
convenient to represent the contribution of the apparent turbulent stress, t^, 
to the total shear stress, t, by an effective turbulent viscosity, and 
similarly the turbulent contribution to the total heat flux, Q, is represented 
conveniently by an effective turbulent conductivity, k<p, such that 
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The turbulent conductivity, kj, is related to the turbulent viscosity, v^, by 
a turbulent Prandtl number, Pr,p, defined as 


Pr^ = p Cp 


(A- 8 ) 


and the boundary layer momentum and energy equations, Eqs. (A-l) and (A-2), 
become 
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In deriving Eq. (A -10) use has been made of the definition of total temperature 
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With the flow laminar, Eqs. (A~3), (A~9), and (A-IO) are solved with v>p = 0 to 
determine the boundary lajrer development. If the flow is transitional or 
turbulent, it is necessary to model Vip and Prtj,. The specification of the 
turbulent viscosity, Vip, and the turbulent Prandtl number, Er-j, is carried out 
through the turbulence model described in detail in the subsequent section. 

After specification of the turbulent viscosity and Prandtl number, 

Eqs. (A“3), (A-9), and (A-IO) are solved by first eliminating from the 
momentum and energy equations by application of the continuity equation. When 
the strearawise static pressure distribution, P(x), is specified, the momentum 
and energy equations, in conjunction with an equation of state and equations 
governing vip and Pr^, form a closed set of nonlinear, parabolic, partial 
differential equations, which can be solved upon specification of boundary 
conditions. The wall and free-stream boundary conditions employed in the 
solution are given by: 

at the wall y = 0 

_ = 0 - 

pv ^ (pv)w, U = 0, T • Tw or : 0 (A-12) 
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at the free-stream, y-»oo 
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Details of the solution procedure are presented in ref. 15. 


The Turbulence Model 

Fully -developed turbulence model. - The fully-developed turbulence model 
originally presented by McDonald and Camarata (ref. 11) for two-dimensional 
incompressible flow, forms the basis for the transitional turbulence model used 
in the McDonald-Fish transitional boundary layer theory and, therefore, at this 
juncture it is useful to describe the model in some detail. The fully-developed 
turbulence model, which is described in greater detail in ref. l4, is based 


) 
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upon a solution of the turbulence kinetic energy equation. The turbulence 
kinetic energy equation is a conservation equation derived from the Navier- 
Stokes equations by writing the instantaneous quantities as a sum of mean and 
fluctuating parts. Hie Navier-Stokes momentum conservation equation 
(i = 1, 2, 3» referring to the three coordinate directions) is multiplied by 
the i'th component of fluctuating velocity and the average of the resulting 
three equations is taken. The three averaged equations are s umm ed to obtain 
the turbulence kinetic energy equation. The derivation of the turbulence 
kinetic energy equation has been given by Favre (ref. 35) for compressible 
flow and approximated by Bradshaw and Ferris (ref. 10 ) to boundary layer flows; 
a derivation and discussion of the turbulence kinetic energy equation used in 
the present theory is given in ref. 15. 


As shown in ref. 15, the boundary layer approximation to the turbulence 
kinetic energy equation is given by 


^ (i pOq* ) pv q*)= - /) W 

advection 


du 
dy 

production 


- *5y (PV T (pv)’ q2) - 
difftMion dissipation 


(A-lU) 


-,...2 


, 1 ^ p 


• P (U* - V’ ) ^ 

Domal atresi production pressure-dilitation 

All calculations reported in the investigation were made with the usual 
assumption of zero pressure-dilitation contribution to the energy balance 
(ref. 10). The turbulence model is developed by integrating Eq. (A-lU) with 
respect to y between the limits y - 0 and y = 6 which leads to 


? 57 


d - - ~2 
p u 


- r 

d* L 


/.# - 

cfi • j - p u' v' 
'o 

du 


du - 

TT Jo P 


p (u'2 . V'*) di P‘ c'y ♦ E 


(A-15) 


where 


E = [*^ q^ (pu -pv) - p'v' (pv) q2] 


(A -16) 
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Following Townsend (ref. 36) and Bradshaw and Ferris (ref. lO) structural 
coefficients a^ and L are introduced, together with a mixing length Z ; these 
scales are defined as 


T" . — 5 2 — y |2 F 

- u V s 0, q' , u' ' Qg , v = Oj q 


w' = ( I - Oj- Oj) q* 


— — 3/2 r '/2 g 

f»(-uv) /L, (-UV) = J ■ 


dy 


(A -17) 


For fully-developed turbulence the structural coefficients a^^, &p, and are 
assumed 'constant having values 0.15, 0.50, and 0.20, respectively (refs. 11 and 
37). As Is discussed subsequently, aj^ becomes a variable in the transitional 
regime. Using Eq. (A -17), Eq. (A-l6) is put in the form 



(A -18) 


where 




(A-19) 



(A -20) 




dr} ) Ug 


du 

dx 


drj 


(A -21) 


where T] is a nondimensional transverse distance y/6^, 6'*’ is an arbitrary 
reference length, and 6 the boundary layer thickness. 
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The left-hand side of Eq. (A-l8) represents the streamwise rate of 
change of turbulence kinetic energy and is derived directly from the turbulence 
kinetic energy advection term. The term PeU^02 represents the integral of 
turbulence production minus dissipation and the normal stress 

production. The terms designated by E are turbulent source terms resulting 
from disturbances imposed upon the boundary layer by the free stream. As 
shown in Eq. (A-l6), E is the sum of two major contributions the first, 
(q2/2)(‘^36/dx-'^) representing the free-stream velocity disturbance (i.e., 
free - s tream t urbulence entrained by the boundary layer) and the second, 
-P’v'+(pv)' ^/2, representing the direct absorption of acoustic energy. In 
subsonic and moderately supersonic flows the acoustic source term is 
negligible compared to the free-stream turbulent velocity disturbance. At 
higher Mach numbers the acoustic term may become the dominant source and, as 
indicated by Pate and Schueler (ref. 3), the acoustic disturbance is the 
dominating disturbance in high Mach number wind tunnels having turbulent 
tunnel wall boundary layers. In the present investigation of the NASA pilot 
model quiet wind tunnel it is the tunnel wall boundary layer itself which is 
under investigation and when this is laminar there is not expected to be an 
acoustic disturbance source. Thus, in the present calculation for the tunnel 
wall boundary layers the acoustic disturbance has been assumed negligible. 

For fully-developed turbulent flow, as in ref. 11, L and I are given by 

= 0. Monh [»fy/(0.l 8)] (A-22) 


~ tonh [ry/l»] (A -23) 


where JL^ is the "wake" value of the mixing length at any particular streamwise 
station. Although Eqs. (A-22) and (A -23) give accurate representations of 
and L through most of the turbulent boundary layer, it is well-known that they 
overestimate the length scales within the viscous sublayer and are somewhat 
inaccurate at low Reynolds numbers. Following McDonald and Fish (ref. l4) the 
experimentally observed damping effect in the viscous sublayer is modeled by 
assuming intermittent turbulence within the sublayer leading to the relation 


32 




-U V 


= r(-u'v’)T = r(iau/<jy4 


(A -24) 


In Eq. (A-24) T is the intermittency factor, ^the dan^jing factor, and the 
subscript T indicates the value with turbulent flow. Obviously, ^ is equal to 
the square root of r . As in ref. l4, the present investigat ion a ssumes that 
the damping distributes normally about mean height y'*’(y''' = y^r/p/v) with a 
standard deviation a leading to the equation 

(A-25) 


where P is the normal probability function; y^ is taken as 23, and a as 8. A 
detailed discussion of the sublayer danping treatment is presented in ref. l4. 

In the present calculations the von Karman constant K was taken to be 0.43. 

In regard to the low Reynolds number effects. Coles (ref. 38) has 
observed and correlated the departure of the mean velocity profile of a flat 
plate turbulent boundary layer from the usual similarity laws known to hold at 
higher Reynolds numbers . Using Coles ' correlation of the mean velocity profile 
in the lew Reynolds number regime, McDonald (ref. 39) integrated the boundary 
layer equations of mean motion to obtain local distributions of turbulent shear 
stress and evaluated the local mixing length distributions from the assumed mean 
velocity distribution and the computed shear stress distributions. Based upon 
these calculations, a low Reynolds number correction for the dissipation length 
of the form 


L = L-« • **P ( - 1 63 In 




9 7 


'] 


(A-26) 


was derived where L<jo/6 is given by Eq. (A -22) . In the calculations presented 
in the present report the dissipation length used was obtained by multiplying 
Eq. (A-26) by the sublayer damping factor 

When numerical values of the structural coefficients an are specified, 

Eqs. (A-22), (A-23), (A-25), and (A-26) are used to represent L and X, and the 
pressure dilitation is either neglected or modeled, the turbulence kinetic 
energy equation, Eq. (A-l8), becomes eui ordinary differential equation with 
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the dependent parameter ioo(x) which is solved in conjunction with the 
boundary layer momentum and energy equations to predict the development of both 
the mean flow field and the turbulent shear stress. 

In addition to including the turbulence kinetic energy equation in the set 
of equations governing the boundary layer develop ment it is necessary to 
specify the turbulent heat flux contribution, -pCpv’T'. As previously 
stated, in the present procedure, v'T' is specified by assuming a turbule nt . 
Prandtl number, Prij, which relates the velocity-temperature correlation, v’T’, 
to the Reynolds stress, u'v', through Eq. (A - 8 ). The turbulent Prandtl number 
distribution used in the present procedure varies with distance from the wall 
in the manner suggested by Meier and Rotta (ref. Uo). At this juncture it 
shou ld be pointed out that an alternative procedure can be used to determine 
v'T’, bas ed upon an easily derive d co nservation equation for either the 
quantity T'2 or the correlation, v’T’, which is similar in form to the 
turbulence kinetic energy equation, Eq. (A-lU). However, to solve this new 
conservation equation it is necessary to assume a universal structure relating 
quantities analogous to dissipation, production, etc. While sufficient 
experimental data exists to allow valid modeling of the required terms for the 
turbulence kinetic energy equation, the existing data does not indicate how 
proper modeling could be carried out for the v'T’ conservation equation. Thus, 
at least for the present, the approach based upon a t urbulent Prandtl number 
appears preferable to an approach based upon the v’T’ conservation equation. 

Transitional turbulence model. - The McDonald -Fish transitional turbulence 
model is based upon a solution of the turbulence kinetic energy equation with 
the structural coefficients modified from their fully-developed turbulence 
values. The following discussion of the model condenses the presentation in 
ref. lU. For fully-developed turbulent flows the structural coefficients a]_, 
& 2 , and a 3 are assumed constant and are set equal to 0.15, O. 50 , and 0.20, 
respectively. Although it is probable that all coefficients vary in the 
transitional regime only the coefficient ai contributes significantly to the 
energy balance when the boundary layer is far removed from separation, as shown 
by Bradshaw (ref. 37). Therefore, the transition model assumes a£ and a^ are 
equal to their fully-developed values and that only a^^ need be modified. The 
McDonald -Fish model modifies a]_ through the introduction of a turbulence 
Reynolds number, R^, defined by 

1/2 

{ - u' vO J 

Rr' — (A-27) 


3 ^ 


Noting that 


(-u'v')* 


2 


i 


( ii )2 , ^ M. 

^ di ‘ (jy 


(A - 28 ) 


the turbulence Reynolds number becomes 


R 


T 



(A -29) 


To be consistent with the integral turbulence kinetic energy equation, 

Eq. (A-18), a layer -averaged turbulence Reynolds number, R^, is introduced as 



(A - 30 ) 


where 6 g, the sublayer thickness, is defined as the location at which the 
laminar stress has fallen to four percent of the total stress (the four percent 
definition gave a sublayer mean temperature in very good agreement with the 
so-called Eckert reference temperature). 

The McDonald -Fish model assumes that the turbulence Reynolds number, R^, 
is the sole variable influencing the development of ap and a relationship 
between a^ and Rj is obtained by considering the development of an incompress- 
ible constant pressure flat plate equilibrium turbulent boundary layer. It 
should be noted that under the assumption that a^ is solely dependent upon R,^ 
it is only necessary to derive a relation for one set of flow conditions to 
obtain a universally valid relationship. Based upon experimental observation, 
it is readily ascertained that for the incompressible, constant pressure 
equilibrium turbulent boundary layer, 0^, E, and h(t>^/hx are negligible which 
leads to the reduced turbulence kinetic energy equation 

din 0 | _ din 8 ^2^ (A- 31 ) 

dR^ ■ dR0 R^ S 


In Eq. (A- 31 ), Cf is the skin friction coefficient and the independent variable 
has been changed from streamwise distance, x, to momentum thickness Reynolds 
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number, R0, using the momentum integral equation. Equation (A-31) has the 
general form of a Bernoulli equation 


din a 
dR® 


- f (R^) + 0| g (R^) 


(A -32) 


which has the general solution 


0 ,: (exp/f dR^)/(C-/gexp/fdRg dR^) (a-33) 


When the mutual dependence between f, g, and R0 is determined for the flat 
plate equilibrium boundary layer, Eq. (A-33) provides the required general 
relationship between ai and Rj. The quantities f and g are evaluated by first 
noting that for the flat plate equilibrium boundary layer 6/6 varies very 
slowly with R0. Neglecting this variation leads to the result 



(A-34) 


Ihe quantity g which is equal to the factor -U026/(Cf0^6) is evaluated 
numerically by Integrating the profiles of Maise and McDonald (ref. 4l). Over 
a wide range of Reynolds numbers the grouping is found to be sensibly constant 
with a value approximately 6.66. Thus, Eq. (A-33) becomes 


0, t Qo (R^/R^^ )/[!•► 6.666 Oo ^R^/Rflo' ')] 


(A -35) 


where the arbitrary constant aQ is the value of ai^ when R0 is equal to ^0.. 

It should be pointed out that at large values of R0/R0Q, a^ asymptotes to the 
fully-developed value of 0.15. The independent variable of Eqs. (A-35) and 
(A- 26 ) is changed from R0 to Rj by using the profiles of Maise and McDonald 
(ref. 4l) which integrate to give 

Rg = 68.1 R,. + 614.3 R^ > 40 (A- 36 ) 
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At low Reynolds numbers good results are obtained using the equation 

_ 0 22 » . , , 

RgrioORr Rt-' (A-37) 


In the intermediate range 1<S^<40, the two distributions, Eqs. (A -36) and 
(A“37)> are joined by a cubic constructed to match the value and slope at the 
join points. Finally, the constant of integration, aQ, is determined on the 
basis of comparison between experiment and theory. Best agreement between 
theory and experiment for low Jfech number, adiabatic wall boundary layers is 
obtained by setting ao equal to 0.0115 when is equal to unity. During the 
assessment of theory as applied to low hypersonic Mach number boundary layers 
Shamroth and McDonald (ref. 15) showed that good results for transitional 
boundary layers could be obtained over a range of Mach numbers and wall to 
free stream temperature ratios by allowing ao to be a function of the wall^to- 
free-stream static temperature ratio, Tw/Tg. The functional relationship was 
chosen so that for incompressible flow the theory would agree with the 
incompressible data of Zysina-Molozhen and Kuznetsova (ref. 42) in the range 
0.5 < Ty,/Tg < 2 .8. For T^j/Tg>2.8 the transition Reynolds number is assumed to 
be independent of wall temperature ratio. The variation of ao with T„/Tg is 
shown in Fig. 4l. The functional relationship based upon incompressible data 
shown in Fig. 4l leads to good comparisons between theoretical predictions and 
experimental data for a variety of boundary layers in the supersonic and low 
hypersonic Mach number regimes (ref. 15). 
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TABLE - NASA SUPPLIED PRESSURE DISTRIBUTION IN 
THE MACH 5 LAMINAR FLOW SLOTTED NOZZLE 


Dimensionless distance 

Dimensionless Radius, 

from throat, 
X = x/r^ 



Approach section and outer wall of slot 


-21.878 

11.5000 

-21.1+00 

11.1+890 

-20.900 

11.4686 

-20.U00 

11.4390 

-19.900 

11.3990 

-19.1+00 

11.3530 

-18.900 

11.3040 

-18.1+00 

11.2530 

-17.900 

11.2020 

- 17 . 1+00 

11.1460 

-16.900 

11.0890 

-16.1+00 

11.0300 

-15.900 

10.9690 

-15.1+00 

10.9080 

-II+.900 

10.8480 

-13.900 

10.7240 

-13.1+00 

10.6550 

-12.900 

10.5816 

- 12 . 1+00 

10.4960 

-11.900 

10.3980 

- 11 . 1+00 

10.2760 

-10.900 

10.1090 

- 10 . 1+00 

9.8810 

- 9.900 

9.6320 

- 9.1+00 

9.3380 

- 6.900 

8.9930 

- 8 . 1+00 

8.6190 

- 3.378 

8.6030 

- 7.878 

8.2160 

- 7.378 

7.8260 

- 6.878 

7.4210 

- 6.378 

7.4210 


Pressure 

Ratio 

P/Po 


1 .0000 


0.9999 

0.9999 

0.9999 

0.9998 

0.9998 





TABLE I. - NASA SUPPLIED PRESSUPE DISTRIBUTION IN 

THE MACH 5 LAMINAR FLOW SLOTTED NOZZLE - Continued 


- 5.878 

6.1360 

0.9997 

- 5.378 

6.1360 

0.9997 

- U.878 

5 .6870 

0.9996 

- 4.378 

5 .2320 

0.9994 

- 3.878 

4.7680 

0.9991 

- 3.378 

4 .2960 

O.99OO 

- 2.878 

3 .8190 

0.9981 

- 2.378 

3.3420 

0.9967 

- 1.878 

2 .8480 

0.9947 

- 1.678 

2.6450 

0.9931 

- 1.638 

2 .6020 

0.9926 

- 1.520 

2 .4730 

0.9910 

- 1.360 

2.2820 

0.9880 

- 1.190 

2 .0600 

0.9828 

- 1.072 

1.8980 

0.9746 

- 0.982 

1.7760 

0.9649 

- 0.908 

1.6780 

0.9536 

- 0.848 

1.6060 

0.9410 

- 0.788 

1.5380 

0.9272 

- 0.720 

1.4680 

0.8485 

- 0.688 

1.4400 

0.7853 

- 0.646 

1 . 4 l 4 o 

0.6968 

- 0 . 6 o 4 

1.4020 

0.6376 

- 0.558 

1.4000 

0.5283 

- 0.504 

1.4100 

0.3343 

- 0.448 

1.4410 

0.1619 

Inner wall of slot 

- 0.900 

1.3740 

0.92716 

- 0.854 

1.3560 

0.88677 

- 0.796 

1.3320 

0.84846 

- 0.736 

1.3100 

0.78525 

- 0.682 

1.2940 

0.69678 

- 0.624 

1 .2840 

0.63764 

- 0.560 

1.2800 

0.52828 

- 0.484 

1 .2800 

0.33427 

- 0.390 

1 .2800 

0.16191 










TABLE I. - NASA SUPPLIED PRESSURE DISTRIBUTION IN 

THE MACH 5 LAMINAR FLCW SLOTTED NOZZLE - Continued 


Transonic throat section 

- 0.9000 

1.374 

0.92716 

- 0.8552 

1.334 

0.91208 

- 0.79^ 

1.2902 

0.89588 

- 0.7U50 

1.253 

0.87863 

- 0.6952 

1.220 

0.86039 

- 0.6U82 

1.1906 

0.84123 

- 0.6034 

1.1650 

0.82120 

- 0.5604 

1.1424 

0.80037 

- 0.5194 

1.1224 

0.77880 

- 0.4796 

1.1044 

0.75654 

- 0.4412 

1 .0884 

0.73365 

- 0.4040 

1.0742 

0.71014 

- 0.3678 

1.0616 

0.68622 

- 0.3324 

1.0506 

0.66179 

- 0.2980 

1.0408 

0.63695 

- 0.2642 

1.0322 

0.61178 

- 0.2312 

1.0246 

0.58633 

- 0.1988 

1.0184 

0.56066 

- 0.1670 

1.0130 

0.53483 

- 0.1356 

1.0086 

0.50892 

- 0.1048 

1.0052 

0.48298 

- 0.0742 

1.0026 

0.45709 

- 0.0442 

1.0010 

0.43130 

- 0.0146 

1.0002 

0.40570 

0.0000 

1.0000 

0.39320 

0.0148 

1.0002 

0.38036 

0.0438 

1.0010 

0.35535 

0.0724 

1.0026 

0.33074 

0.1010 

1.0050 

0.30662 

0.1280 

1.0084 

0.28909 

0.1572 

1.0124 

0.26015 

0.1850 

1.0174 

0.23839 

0.2400 

1.0300 

0.19603 

0.2944 

1.0460 

0.15789 





TABLE I. - NASA SUPPLIED PEESSUEE DISTRIBUTION IN 

THE MACH 5 LAMINAR FLOW SLOTTED NCEZLE - Continued 


Supersonic section 

0.3512 

1.0654 

0,15700 

0.3992 

1.0824 

0.14861 

O.U558 

1.1036 

0 .l 40 l 8 

0.5218 

1.1286 

0.13162 

0.5972 

1.1438 

0.12287 

0.6826 

1.1924 

0.11437 

0.7766 

1.2306 

0.10628 

0.8788 

1.2728 

0.09851 

0.9902 

1.3190 

0.09087 

1.1118 

1.3696 

0.08355 

1.2420 

1.4238 

0.07667 

1.3804 

1.4812 

0.07027 

1.5256 

1.5410 

0.06448 

1.6778 

1.6028 

0.05916 

1.8358 

1 .6662 

0.05435 

1.9992 

1.7310 

0.05000 

2.3440 

1.8640 

0.04257 

2.7080 

1.9996 

0.03622 

3.0940 

2.1368 

0.03091 

3.5020 

2.2740 

0.02660 

3.9300 

2.4106 

0.02302 

4.3760 

2.5464 

0.02006 

4.8400 

2.6804 

0.01758 

5 .3240 

2.8122 

0.01549 

5.8280 

2.9418 

0.01373 

6.3520 

3 .0648 

0.01222 

6.3960 

3.1914 

0.01093 

7.4600 

3.3112 

0.00982 

8.0380 

3.4280 

0.00887 

8.6360 

3.5408 

0.00805 

9.2580 

3.6488 

0.00732 

9.8900 

3.7542 

0.00670 

10.5440 

3 .8548 

0.00614 

11.2220 

3 .9500 

0.00566 

11.9040 

4 .0436 

0.00523 

12.6200 

4.1294 

0.00484 

13 .3360 

4. 2140 

0.00451 

14.0900 

4.2900 

0.00420 
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14.8380 

15.6240 

16.4060 

17.2240 

18.0380 

18.8920 

19.7440 

20.6040 
21.4860 
22.3840 
23 .2900 
24.2160 
25.1500 
26.0920 
27.0500 
28.0140 
28.9860 
29.9640 
30.9500 


NASA SUPPLIED PRESSURE DISTRIBUTION IN 
THE MACH 5 LAMINAR FLOW SLOTTED NOZZLE 


4.3664 

4.4334 

4.5005 

4.5582 

4.6158 

4.6634 

4.7106 

4.7550 

4.7928 

4.8258 

4.8564 

4.8806 

4.9026 

4.9214 

4.9354 

4.9475 

4.9568 

4.9636 

4.9680 


Concluded 


0.00394 

0.00369 

0.00348 

0.00328 

0.00311 

0.00295 

0.00281 

0.00269 

0.00257 

0.00247 

0.00237 

0.00229 

0.00222 

0.00215 

0.00210 

0.00205 

0.00200 

0,00197 

0.00194 
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TABLE II. - NASA SUPPLIED PRESSURE DISTRIBUTION FOR 
INLET DUCT AND APH^OACH SECTION BOUNDARY 
LAYER CALCULATION 


Dimensionless Distance 
^ From Throat , X = 

Contour Distance From 
Last Screen, S 

m. ft. 

Dimensionless 

Radius, 

Pressure 

Ratio 

-165.878 

0.0 

0.0 

11.5000 

0.9999770 

- 21.878 

1.3288 

6 .000000 

11 .5000 

0.9999770 

- 17.400 

1.8345 

6.018718 

11.1460 

0.9999673 

- 15.900 

1.9050 

6.250103 

10.9690 

0.9999714 

- 1U.900 

1.9178 

6.292076 

10.8480 

0.9999583 

- 13.900 

1.9306 

6.334054 

10.7240 

0.9999432 

- 12.900 

1.9434 

6.376114 

10.5810 

0.9999342 

- 11.900 

1.9563 

6.418391 

10.3980 

0.9999304 

- 11.400 

1.9694 

6.461433 

10.2760 

0.9999281 

- 10.900 

1.9763 

6.483864 

10.1090 

0.9999240 

- 10.400 

1.9833 

6.506949 

9.8810 

0.9999169 

- 9.900 

1.9905 

6.530670 

9.6320 

0.9999052 

- 9.400 

1.9981 

6.55541 

9.3380 

0.9998876 

- 8.900 

2.0059 

6.581074 

8.9930 

0.9998619 

- 8.400 

2.0063 

6.582213 

8.6190 

0.9998309 

- 8.378 

2.0142 

6.608265 

8.6030 

0.9998293 

- 7.878 

2.0222 

6.634649 

8.2160 

0.9997898 

- 7.378 

2.0326 

6.6688233 

7.8260 

0.9997433 

- 8.900 

2.0977 

6.688233 

7.4210 

0.9996898 

- 6.378 . 

2.0469 

6.715564 

7 . 004 o 

0.9996292 

- 5.878 

2.0553 

6.743150 

6.5720 

0.9995614 

- 4.878 

2.0638 

6.770971 

6.1360 

5.6870 

0.9994841 

0.9993846 

- 4.378 

2.0723 

6.799056 

5 .2320 

0.9992465 

- 3.878 

2.0810 

6.827351 

4.7680 

0.9990510 

- 3.378 

2.0897 

6.355887 

4.2960 

0.9987087 

- 2.878 

2.0984 

6.884608 

3.8190 

0.9980443 

- 2.378 

2.1072 

6.913401 

3.3420 

0.9968761 

- 1.878 

2.1160 

6.942441 

2.8480 

0.9946463 

- 1.678 

2.1196 

6.954235 

2.6450 

0.9930038 

- 1.638 

2.1204 

6.956646 

2.6020 

0.9925850 

- 1.520 

2.1226 

6.963898 

2.4730 

0.9910844 

- 1.360 

2.1257 

6.973979 

2.2820 

O.9S8197O 

- 1.300 

2.1269 

6.977909 

2.2036 

0.9867291 








TABLE II. - NASA SUPPLIED PRESSURE DISTRIBUTION FOR 
INLET DUCT AND APPROACH SECTION BOUNDARY 
LAYER CALCULATION - Concluded 


- 1.190 


6.985453 

2.0600 

0.9827085 

- 1.120 


6.990462 

1.9638 

0.9785564 

- 1 . 0?2 

2,1317 

6.993906 

j. »898 o 

0.9747525 

- 1.030 

2.1326 

6.996885 

1.8410 

0.9652792 

- 0.982 

2.1337 

7.000220 

1.7760 

0.9652792 

- 0.940 

2.1346 

7.003150 

1.7204 

0.9595346 

- 0.908 

.2.1352 

7.005320 

1.6780 

0.9542777 

- 0.870 

2.1360 

7.007820 

1.6324 

0.9465237 

- 0.848 

2.1364 

7.009240 

1.606 

O.94Q6963 

- O.820 

2.1369 

7.611020 

1.5742 

Oi^o 6 ce 7 

- 0.788 

2*1376 

7.013090 

1.52^ 

0.9076510 

- 0 . 7 C'' 

2.1377 

7.013590 

1.5298 

0.9076510 

- 0.758 

2.1381 

7.014920 

1.5670 

0.8888883 

- 0.740 

2.1396 

7.01980 

1.4886 

0.8683142 

- 0.720 

2.1388 

7.017170 

1.4680 

0.8393871 

- 0.710 

2.1390 

7.017740 

1.4592 

0.823443 

- 0.688 

2.1394 

7.018890 

1.4400 

0.7873495 

- 0.666 

2.1397 

7.020020 1 

1.4264 

0.7470208 

- 0.646 

2.1400 

7.021010 ^ 

1 . 4 l 40 

0.7079584 

- 0.630 

2.1402 

7.021760 

1.4094 

0.6761727 

- 0.604 

2.1406 

7.022910 

1.4020 

0.6239762 

- 0.585 

2.1408 

7.023710 

1.4012 

0.5853385 

- 0.568 

2.1410 

7.024420 

1.4004 

0.5495445 

- 0.558 

2.1411 

7.024840 

1.4000 

0.5277362 

- 0.552 

2.1412 

7.025090 

1.4012 

0.5145390 
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TRANSITION REYNOLDS NUMBER 


u*fn 

THEORY OF vin DRIEST ft BLUMER.“JJ- » 0.01'IREF. 18) 

UARL PREDICTION, — ^=0.012 
DATA OF FEINDT, -^=0.012 iREF. 17) 



Figure 1. — Comparison between measured and predicted transition Reynolds number. 
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DIMENSIONLESS SHAPE FACTOR, H MOMENTUM THICKNESS REYNOLDS 

MIXING LENGTH -^„/5 NUMBER. Ra x 



















MENSIONLESS SHAPE FACTOR - H MOMENTUM THICKNESS REYNOLDS 



















CURVES AND EXPERIMENTAL DATA FROM JONES AND LAUNDER (REF. 22) 

O- PREDICTIONS UNDER THE CONDITION OF JONES AND LAUNDER, -S2. = o.Ol 

U 

A-NEW DATA OF JONES AND LAUNDER (REF.22) 

TURBULENT THEORY 


LAMINAR THEORY 




Figure 9. — Comparison between predictions and measurements of momentum thickness 
Reynolds number and shape factor for sink flow boundary layers. 
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STANTON NUMBER 


Q DATA OF MORETTI AND KAYS RUN 12 (REF. 29) 
PREDICTION ( — = 0.01) 

PREDICTION FULLY TURBULENT 




DIMENSIONLESS STREAMWISE DISTANCE, x = X/l^ef 


Figure 12. — Comparison between predicted and experimental results for heat transfer in an accelerated flow. 


WALL 


DIMENSIONLESS MIXING LENGTH, £oe/« STANTON NUMBER 


o 


DATA OF MORETTI AND KAYS RUN 41 (REF. 29) 





Figure 13. — Comparison between predicted and experimental results for heat transfer in an accelerated flow. 






DIMENSIONLESS MIXING LENGTH.^oo/® STANTON NUMBER 


Q DATA OF MORETTI AND KAYS RUN F-2 (REF, 29 
PREDICTION ( ^ - 0.01) 

— PREDICTION FULLY TURBULENT 




Figurs 14. — Comparison between predicted and experimental results for heat transfer in an accelerated flow. 
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<TwALL"^oo* 



DIMENSIONLESS NOZZLE RADIUS, r/r, 



DIMENSIONLESS DISTANCE FROM THROAT. X= x/r^ 


Figure 15. — Gaiculated streamlines md sonic line for the transonic section of the NASA Mach 5 nozzle 
with slot fared in— throat and supersonic section follow NASA contour. 


63 




-3.0 -2.0 -1.0 0 1.0 2.0 3.0 

DIMENSIONLESS DISTANCE FROM THROAT, X= x/r^ 


Figure 16. - Calculated streamlines and sonic line for the transonic section of the NASA Mach 5 nozzle 
using estimated slot stagnation streamline contour— throat and supersonic section follow 
NASA contour. 
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DIMENSIONLESS NOZZLE RADIUS, r/r. 


CONTOUR A 


.CONTOUR B 



DIMENSIONLESS DISTANCE FROM THROAT.X= x/r^ 


Figure 17. — Comparison of calculated sonic lines in the NASA Mach 5 nozzle 
using different approach section contours. 



PRESSURE RATIO. P/P, 


----—CONTOUR A 
CONTOUR B 

— — — — — —NASA DISTRIBUTION 



Figure 18. - Comparison between calculated and NASA wall pressure distribution in the transonic 
section of the NASA Mach 5 nozzle — calculations based on two different 
approach section contours. 
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PREDICTION USING UAC TRANSONIC FLOW CODE 
NASA SUPPLIED PRESSURE DISTRIBUTION 
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throat of the NASA Mach 5 nozzle. 



















CALCULATION INITIATED AT SLOT - THROAT AT X* 0 


u"™ 

!2«0.01 

U) 

JLin. -0.06 

Ui 


■gj— - 0.03 



Figure 25. — Acceleration paranwter and predicted incompressible shape factor in the transonic and 
supersonic sections of the NASA Mach 5 laminar flow slotted nozzle. 
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SHAPE FACTOR, H MOMENTUM THICKNESS , 







BOUNDARY LAYER THICKNESS, 5 DIMENSIONLESS MIXING LENGTH, C „/5 
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Figure 27. - Predicted dimensionless mixing length and boundary layer thickness in the transonic 
and supersonic sections of the NASA Mach 5 laminar flow slotted nozzle. 


75 



INCOMPRESSIBLE SHAPE FACTOR, H* ACCELERATION PARAMETER, 


CALCULATION INITIATED AT SLOT-THROAT AT X =0 



DIMENSIONLESS DISTANCE ALONG NOZZLE AXIS, x = X/L,-ef 


Figure 28. — Effect of heat transfer on acceleration parameter and incompressible shape factor in the 
transonic and supersonic sections of the NASA Mach 5 laminar flow slotted nozzle. 
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SHAPE FACTOR, H MOMENTUM THICKN ESS, 9 








BOUNDARY LAYER THICKNESS, 6 DIMENSIONLESS MIXING LENGTH 










AT SLOT-THROAT AT X * 0 


2- = 0.03 


FLOW CONDITIONS 
pQ = 10.2 atm (150 psia) 

Tq =65.5°G(150°F) 

Uj = 1 1 76 cm/sec (40 fps) 

— '•’wall ' ADIABATIC 

■'■yyj^LL =-17.8°C (0°F) 

Tyu.i L =-73.3°C(-100°F) 




SHAPE FACTOR, H MOMENTUM THICKNESS, 0 


CALCULATION INITIATED AT SLOT-THROAT AT X = 0 



Figure 32.- Effect of heat transfer on momentum thickness and shape factor in the transonic and 
supersonic sections of the NASA Mach 5 laminar flow slotted nozzle. 







DIMENSIONLESS DISTANCE ALONG NOZZLE AXIS, x = X/L^gf 


Figure 33. — Effect of heat transfer on mixing length and boundary layer thickness in the transonic 
and supersonic sections of the NASA Mach 5 laminar flow slotted nozzle. 
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INCOMPRESSIBLE SHAPE FACTOR, H‘ ACCELERATION PARAMETER 


CALCULATION INITIATED AT SLOT - THROAT AT X = 0 



Figure 34. — Effect of stagnation pressure on acceleration parameter and incompressible shape factor 
in the transonic and supersonic sections of the NASA Mach 5 laminar flow slotted nozzle. 
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SHAPE FACTOR, H MOMENTUM THICKIMESS, ^ 
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BOUNDARY LAYER THICKNESS, 5 DIMENSIONLESS MIXING LENGTH, 





MOMENTUM INCOMPRESSIBLE 

THICKNESS, « SHAPE FACTOR, H* ACCELERATION PARAMETER, K ="i^ 


[IW*] 


CO 
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I 8 





DIMENSIONLESS DISTANCE ALONG NOZZLE AXIS, X » X/L^gf 


Figure 37. — Effect of nozzle geometry on acceleretion parameter, incompressible shape factor and 
momentum thickness in the transonic and supersonic sections of the NASA Mach 5 
laminar flow slotted nozzle. 
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BOUNDARY DIMENSIONLESS SHAPE FACTOR, 

LAYER THICKNESS, 8 MIXING LENGTH, 


FLOW CONDITIONS 
T(,-65.5°C(160°F) 

Pq = 10.2 atm (150 psta) 
ADIABATIC 

U j * 11 76 cm/sec (40 fps) 

NASA CONTOUR 

SHARP CORNER 

THROAT 


DIMENSIONLESS DISTANCE ALONG NOZZLE AXIS, x = X/L^ef 


Figure 38. - Effect of nozzle geometry on shape factor, mixing length and boundary layer thickness 
in the transonic and supersonic sections of the NASA Mach 5 laminar flow slotted nozzle. 





MOMENTUM THICKNESS, 6 INCOMPRESSIBLE SHAPE FACTOR. ACCELERATION PARAMETER, 


CALCULATION INITIATED AT THROAT 
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DIMENSIONLESS DISTANCE ALONG NOZZLE AXIS, x = X/L^gf 


Figure 39. - Effect of Mach number on boundary layer development in supersonic nozzles. 





BOUNDARY LAYER THICKNESS, 6 DIMENSIONLESS SHAPE FACTOR, H 

MIXING LENGTH. 






WALL TEMPERATURE 
FREE-STREAM TEMPERATURE' 


Figure 41. — Variation of structural coefficient with temperature ratio. 
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